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Аннотация 

Введение. Для описания работы температурных пьезокерамических конструкций используется теория 
термоэлектроупругости, в которой математическая модель сформулирована в виде системы несамосопряженных 
дифференциальных уравнений. Сложность ее интегрирования в общем виде приводит к исследованию задач в 
несвязанной постановке. Это не позволяет оценить эффект влияния электроупругих полей на температурное 
поле. В литературе не представлены исследования данных задач в трехмерной связанной постановке, в которых 
были бы построены замкнутые решения. При этом проведение именно таких исследований позволяет понять 
картину взаимодействия механических, тепловых и электрических полей в конструкции. Поэтому целью 
представленного исследования стало построение нового замкнутого решения связанной задачи для 
пьезокерамической круглой жестко закрепленной пластины, позволяющее качественно оценить взаимное 
влияние термоэлектроупругих полей в данной электроупругой системе. 

Материалы и методы. Объектом исследования является пьезокерамическая пластина. Рассматривается случай 
нестационарного изменения температуры на ее верхней лицевой поверхности при учете конвекционного 
теплообмена нижней плоскости с окружающей средой (граничные условия 1 и 3 рода). Индуцируемое в 
результате образования температурных деформаций электрическое поле фиксируется путем подключения 
электродированных поверхностей к измерительному прибору. Задача термоэлектроупругости включает 
уравнения равновесия, электростатики и нестационарное гиперболическое уравнение теплопроводности. Она 
решается обобщенным методом конечного биортогонального преобразования, позволяющего построить 
замкнутое решение несамосопряженной системы уравнений. 

Результаты исследования. Построено новое замкнутое решение связанной осесимметричной задачи 
термоэлектроупругости для круглой пластины, выполненной из пьезокерамического материала, являющееся 
более точным, по сравнению с тем, которое было разработано при решении задач в несвязанной постановке. 
Обсуждение и заключение. Полученное решение начально-краевой задачи позволяет определить 
температурное, электрическое и упругое поля, индуцируемые в пьезокерамическом элементе при произвольном 
температурном осесимметричном внешнем воздействии. Проведенные расчеты дают возможность рассчитать 
размеры сплошных электродов, которые позволяют повысить функциональные возможности пьезокерамических 
преобразователей. Численный анализ результатов помогает выявить новые связи между характером внешнего 
температурного воздействия, процессом деформирования и величиной электрического поля в пьезокерамической 
конструкции. Это дает возможность обосновать рациональную программу экспериментов при их 
проектировании и значительно сократить объем натурных исследований. 


Ключевые слова: задача термоэлектроупругости, связанная осесимметричная задача, жестко закрепленная 
пластина, биортогональные конечные интегральные преобразования 
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Введение. Для повышения функциональных возможностей пьезокерамических датчиков [1-3], основанных 
на взаимном влиянии термоэлектроупругих полей, используются различные математические модели. Для более 
точного учета эффекта связанности данных полей возникает необходимость в построении замкнутых решений. 
Для решения систем исходных несамосопряженных дифференциальных уравнений применяются некоторые 
упрощения. Таким образом, задачи могут быть рассмотрены в несвязанной постановке, либо в задачах 
рассматриваются и анализируются элементы, у которых имеется вырожденная геометрия. В [4, 5] 
рассматривается несвязанная стационарная задача для длинного электроупругого цилиндра, а работа [6] 
посвящены анализу тепловых напряжений в полой сфере. Исследования [7,8] связаны с определение 
температурного поля в пьезокерамической оболочке и круглой пластине при решении несвязанных задач. 
Связанные динамические задачи для однородного пьезокерамического слоя, а также динамические задачи в 
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связанной постановке для градиентно-неоднородного пьезокерамического слоя рассматривались в трудах [9, 10], 
в статьях [11, 12] анализировались поля в неограниченной среде. В работах [13, 14] рассматривается длинный 
полый цилиндр и анализируются термоэлектроупругие поля. 

В настоящее время в литературе не описаны результаты построения замкнутых решений упомянутых 
нестационарных задач в трехмерной связанной постановке. Поэтому в настоящей работе рассмотрена круглая 
пластина, выполненная из пьезокерамического состава и имеющая жесткое закрепление, для которой получено 
новое замкнутое решение задачи термоэлектроупругости. Использование ограничения на скорость изменения 
температуры на ее лицевой поверхности [10] позволяет не включать инерционные характеристики исследуемой 
системы и применять уравнения равновесия в расчетных соотношениях. 

Материалы и методы. В процессе решения использовали обобщенное конечное биортогональное 
преобразование, позволяющее понизить размерность несамосопряженной системы уравнений и, существенно 
упрощая исследования в пространстве изображений, построить замкнутое решение. 


Математическая модель. Рассмотрим некоторую область © : {0 <и, <5,0<0<21,0<2, <й"*}, которую в 
цилиндрической системе координат (7,,0,2.)занимает круглая в плане пьезокерамическая сплошная пластина. 
Для рассматриваемой задачи можно использовать произвольные температурные граничные условия. Однако для 
определенности решения на верхней (2, =0) лицевой поверхности учитывается изменение температуры 
©! (г...) при заданной температуре внешней среды 9” на нижней (2,=Й"”) плоскости (1, — время). 


Цилиндрическая теплоизолированная поверхность жестко закреплена: отсутствует радиальная компонента 
вектора перемещений и угол поворота, а в вертикальной плоскости зафиксирована ее нижняя часть. Нижняя 
Плоскость рассматриваемой круглой пластины имеет заземление. Лицевые электродированные плоскости 
пластины подключены к измерительному прибору. Расчетная схема пластины представлена на рис. 1. 


2» 


Рис. 1. Расчетная схема пластины 


Математическая формулировка рассматриваемой задачи в безразмерной форме для аксиально 
поляризованного пьезокерамического материала с гексагональной кристаллической решеткой состава 6 мм 
имеет вид: 
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а =а-Б/Л, ©"(т.,2.,1,) — приращение температуры в размерной форме; И”(и,, 2.,1„), *(.,2., 1) 


ф*(".,2.,й,) — компоненты вектора перемещений, потенциал электрического поля; 
0: (Г, 2, 1), О» (Г, 2, ) — компоненты тензора механических напряжений; Ш) (г, 2, {) — радиальная компонента 
вектора индукции электрического поля; Л, К, о, — коэффициенты теплопроводности, объемной теплоемкости и 


линейного температурного расширения; фо(",Г) — электрический потенциал, индуцируемый на верхней 


лицевой поверхности; ул, ®; — компоненты тензора температурных напряжений и пирокоэффициентов 


(= 1,3, уз = Слю,); е15, еза, езз, 811, &з —^ пПьезомодули и коэффициенты диэлектрической проницаемости; 
9*=Т-Т. Т, То— текущая температуры и температура первоначального состояния тела; В„е/— время 
1 


релаксации; а — коэффициент теплоотдачи, 9. —_ известная скорость изменения температуры; У = ря + 
гг 


Для определения потенциала электрического поля, индуцируемого в процессе деформирования, на верхней 
лицевой поверхности фу (›, /), в случае подключения электродов к измерительному прибору с большим входным 


сопротивлением, используется дополнительное граничное условие: 


д 
5 [| 2-8 =0 (5) 
(5) 


где р: (г, 2, 0 —_ аксиальная компонента вектора индукции; $ — площадь поверхности. 
Построение общего решения. Для выполнения условия закрепления цилиндрической поверхности пластины 


в вертикальной плоскости вводятся новые функции " (и, 2, 1), И (В: 
И (*, 2,1) = (И-+и(к, 2. 1), (6) 
что позволяет сформировать краевую задачу относительно функций 0, у», ф, ©, которая исследуется методом 


конечных преобразований Фурье-Бесселя: 
1 


Ин (или) = [ (ездил (д) (7) 


0 


[№н (п,2,1),фи (п,2,ё),М№н (п, 2,1} = [>(..24).96.24).0 (0.2.1) Ло (ик) ак, 


н (П,2,1 . 
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где /„ — положительные нули функции .Л| ( т (" = 0,0; Ло = 0}; У, (...) — функции Бесселя. 
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Здесь необходимо отметить, что для удовлетворения последнего граничного условия (2) необходимо принять, 
что первоначальная температура пластины Ту равна температуре внешней среды 9” и функция приращения 
температуры на верхней лицевой поверхности (1, #) =0. Данные допущения без большой погрешности 
позволяют считать, что на цилиндрической поверхности пластины ®(1, 2, # = 0. 

В результате использования алгоритма преобразования в области изображений получается следующая 
начально-краевая задача: 


д?и ди 0ф 
—72ин +а Нл ч+ал И + Мн =0, 9 
Лин та 22? 2] ра 3/ ра ЛУ Н (9) 
д? ди 07ф ОМ 
дан ааа -а5]2фн+ Ш =, 
Л \н т 94 02? 2] п р @5.)]„фи +46 2? @7 ри 
д? ‚ ди - и | [м 
фи аз Фи +9) “-аюЛлин ап и чар), Мн +ав—"=0, 
02 02 02 
92мм д г“ ди 0ф 
:2 Н . Н Н 
Мн + + Мн +а пин + а = 0; 
ЛИН 0? [2 =] Н «[. Н 02 15 ро 
д д д 
2=0,#; аб лин +а4 он + 46 фи а7Мн = 0; сн Лами = 0, (10) 
02 02 02 


ом 
Фн|--о = Фон, Фи =0, Мило = б9н, Св, вом =а Эн; 
[в=Й 


02 
Э\ин, 
о орон, = 0, (11) 
Ге |/=0 
ди | Ом . 
и оН› и = №н; 


= 
О = О и=0 


\ АЙ’, (1 
где бронь Я] [43-8 (0),- | 
(=0 
0 


На следующем этапе решения введение функций Ин(п,2,1), И’н(и,2,6),Фн(п,2,6),Он(и,2,6) при 
использовании следующих соотношений: 


Ин (п, 2,1) =Н, (п,2,6)+Он (п,2,1), Ун (п,2,1) =Н, (п,2,1)+Ин (п,2,1), (12) 


фн (п, 2,1) =Н} (п,2,1)+Фн (п,2,1), Мн (п,2,1) =Н. (п,2,1)+Он (п,24), 
позволяет приведения условий (9) к однородным. 
Здесь {Н/..Н4}={Н!..Н1}+{)(2)..Ло (2)} фона), НЕ.Н} = {9 (2).„Ли (2) ине) +5 (2). (2)} Зы, 
Л (=)...Л› (2) — дважды дифференцируемые функции. 
Подстановка (12) в (9) — (11) при удовлетворении условий: 


1 


‚ © | (ик) ак. 


ОН Н 
2=0,й; ав /,Н, +а4 + @6 а7Н4 =0; о г 
02 02 02 


=0, (13) 


ОН 
Ну,-о = Фон, Нз-ь = 0, Н4--о = бн, [9+ вой, =а Эн, 


[в=й 


позволяет сформулировать следующую задачу: 


0? ‚ 01’ ‚ д 
-ЛИн +а: = Я Ли — чан =Ё, (14) 
0217. 90 0ф 90 
.2 Н . Н ‚2 Н Н 
—а Ли Й’н + а а5]иФн + =Р, 
Ян т а4 2? а2] ра @5/„Фн т @6 02? @7 92 2 
д? ‚ ди 027, д 
фин - а а + @9./„ = ао Ин + ап те + а12 „Он +аз РЗ 
920 а ОИ’ 06 
— 20+ -| —+В— +а4| Ин+- И |-а5- | =Е,, 
Он 922 [2 Ро Он +а4| Он ри о 4 
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Ин 
02 


Йн Офн 


2= 0,1; аб „Он +а4 - +а6 а7Он =0; /.Й’н = 0, (15) 
02 02 


д 
Фи|-=0 < Фи|== — 0, Оно = 0, [нато — 0; 


[ЕЙ 
1=0; Чи = Ион; Ин = Шон; Фн = Фон; Он = Он, (16) 
о =Оон, _ = Рон __ = Фон, ой — Он; 
где 
Е =ЛН, -а, р +а>], та аз» ее На, 
Е, =а,,Н —а4 р а» т +45 „Ну -—@5 р + 47 а , 
В =-Ну+аз и», наи, —ап р —а2 „На -@з р › 
В = На о +2 9]. ани, а о 
Обн = Ни, Тн = "он — Нэу-о› Фон = —Нзу-о› Обн — —Н 4-0» 
Оон = ый 2 Й)н = он [7 ‚ Фон = [о 2 Обн = он и 


Используя биортогональное конечное преобразование (КИП) [15], получаем решение задачи (14) - (16). На 
сегменте [0,И] вводится КИП с неизвестными компонентами вектор-функций преобразований 
Ки, 2). : „КаОм, 2), М! (Ци, 2). : „М 2): 


й 
бб.) = | |, +2 [ни + их | Он | (17) 
7 
0 


02 


{м (иы›2), № (иы›2), М (иы,2), № (ны,2)} 
Кы[ 


т 


{Чн,Тн,Фн,Он} — У би.) 


1 


[к 


Г 
рь | Ки, 2 (инь, 
0 


где Л», И„— собственные значения соответствующих задач относительно компонент вектор-функций 
КИП (% = 1...4). 
В ходе преобразований получаем задачу для определения трансформанты С(и, Аш, 8): 


И. а , — 
раки ббьльи) = ибо (198 п=0,>}, (18) 


= 0, С(^„,п,0) = Он, 


= Со (^„,п), (19) 


решение которых будет иметь следующий ВИД: 


С(п,^„,Ё)= (ты — Тя )' (С — Сотоы )ехр (ть) = (бо — бот ехр (тьы!) + 


1 

+ в [2 (п. т) ехр(тзь (#— т)) - ер(т (#— т) Чт}, (20) 
0 

кроме того, две однородные задачи относительно компонент К1(Аи, 2)...К4(№ 2), 
а?Ки АК»; ак»; 
;2 |т ; 2т . Зт 2 ; 
СЕ „Ки +а Я2 Ли Я. +Льа „К т = 0, 21 
Л" Вл 1 а 2] ря 9] 2. 14/54 (21) 
ь а?К т . АК т ; 4?К т ак т 
—а,.Коы +а4 и +4), — "ар Кзи +ап > Ааа р =0, 
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Ку аКь ЯК АК, 
;2 Зт ; т :2 2т 2 4т 
„К па Яз.Л» а „К т +а +Л ма = 0, 
Л Вз 8 т 3 Е 5/2 6 РЕ. 15 га 
2 
(2.2, РК -е а. + Ки ча — "+ар,Кзы фаз" =0; 
АК; АК з; 
2=0,1, аб. „Ки + а4 — а — Ара 5Каи = 0, Кзы-о = Кзы-ь = 0, (22) 
АК; АК 4; 
ти Кэт ых КЗ 5. 0, К4и=-0 а 0, г + во. = 0; 
42 а 42 =Й 
и М! (Ци 2)...М№(иш, 2): 
;2 а - 1т : а 2т Зт ; 
Ут М: +а: 2 2. п + аз + М т = 0, (23) 
4 42 
а? М»: а? №: АМ 
;2 2т . т :2 Зт 4т 
а Ми а. +ча —а5 „№ +а а = 0, 
12 4 р 2] т 5./п {Уз 6 Е? 7 р 
а? М, АМ; ЯМ. АМ№л, 
;2 Зт : | :;2 2т ; 4т 
М —а + а9] п ао) М чан +а12 „Мы +а = 0, 
Лт Из а 9] 2 10./п 12 я 12.]п1\ 4 13 ра 
а Ма .. 
АМ + в На М4 и аа. Мат + @14 р = @15 зи = 0; 
Е 42 42 
АМ»; АМ; ОМ, 
2 0, й, ЧА + @4 2 + @6 ри а == 0, т Мот — 0, (24) 
42 4 02 
ОМ 
М0 = М щей = 0, М4т-=0 = 0, - + вом.) = 0; 
02 я 
где 
й 
Ен (п, „,Е) — [кк +Е>Ко + ЕК зи + К 4 Е, 
0 
5 Е | р а{Тн,Тн} а {Фон,Фон} 
{би бн} = | {Оон, Он +ац Л» {Оон,Оон} + 415 4т 4, 
и 4 42 


пи, Тэи — корни характеристического уравнения: Вт? +т„ +А2. =0. 

Построенная однородная задача (23), (24) относительно функций М (щи 2)...М№ (и, 2) инвариантна исходным 
расчетным соотношениям (14), (15). 

Системы (21), (23) приводятся к следующим уравнениям относительно К>(\, 2), №( щи, 2): 


а а° а“ Я 
[5 Нет 42 + ет 92“ + ет 42 + @4т (Ки, Мы} =0. (25) 


Коэффициенты е1и...е4ш В статье не приведены в связи ограничением ее объема. 
В уравнении (25) левую часть разложим на коммутативные сомножители, представленные ниже: 


ты а? а“ а? 
2 2 2 2 и 
[© В Ат 2 а Аи Е“ +3 2 +714 т {Ки Мы} =0, (26) 
ел; 
2 и 4 > 
где А, =\/Вы, Аи = Эш, тзи =е +В + Эт, Тди = РА, В „,.5и — Действительные положительные 
1т-2т 


корни следующих характеристических уравнений: 


4 3 2 = 
Ви + Вы ве е›иВи БЕ езиВы а ди = 0, 


ел: 
3 2 2 4т 
Ку — (2 +Вш 52. +(е»Вьь +В: чем бы =. = 0. 
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При исследовании круглой жестко закрепленной пьезокерамической пластины общий интеграл 
уравнений (26) будет иметь следующий вид: 
{Ки Мот } — {Ра ‚Е } ехр (Аш2) + {Р.„ ‚Еэт } ехр (-Аш2) а О ‚Езт } 5т (45и2) зе (27) 


+{Р. ыы С05 (А›и2) + {Р,Еби } 5т( Аз 2) + (Ре, Евы} 05 (Ази2) + 
+{Р,Ет} 5т (442) + {Ри Еви} со5( Аи), 


0,5 0,5 
2 4 2 2 Г 4 2 
— [= +73 —4т4и А — [= —\/Тз —4т4 
› 9 4т ы 


где 


2 


Здесь следует отметить, что условие о действительных положительных значениях коэффициентов 
Вт, Ут, Алт...А4ди выполняется для большинства конструкций, выполненных из пьезокерамического материала. 
В противном случае просто меняется структура формул (26), (27). 

Учитывая, что ранее были получены связи в результате приведения (21), (23) к (25), получаем выражения для 
функций К\(Аш, 2), Кз(№, 2), Кам, 2), Мл, 2), Мм, 2), М4, 2). 

Подстановка К1(Ли, 2)...Ка(№и, 2), М(№, 2)...М(щ, 2) в условия (22), (24) позволяет определить постоянные 
Рут... .Озт, Елт ...Езт И собственные значения Ли, Ши. 

Итоговые выражения функций И (п, 2, 1, И (п,2, 1), ф (п,2, 0), © (п,2,6) получим, применяя формулы 
обращения (17), (8). Тогда, с учетом (6), (12), имеем: 


ме . 
[#8 2,1) =2 > : (7 ин п, 2,4) У С(п, ^ М, ( Ни, 2 Ук» | (28) 
п=1 Ло (1, 1 


т. | (и,21)+ У `б(и Аим, (ии,2) Ки -| 
уй 1 


Ло (у, 
т ( 
ВХР: о пил аольмфыдр -| 


До +6» АМ (и, "| 


п=0 Ло ТР РЕ] 
Функции Л (2)... №2(2) вычисляются из условия упрощения Ё1... Ад при удовлетворении условий (13): 
9°Н дН дН 
;2 1 . 2 й 3 . 
„Ну. -а + а>] п ЯЗ] „На =0, 
ЛЮ 1 1 2? 2] 22 37 22 вЫ 4 
9*Н ОН 9?Н ОН 
а. Н> —а4 > аа Нз-— ав > +а-——=0, 
02 де 
0?Н ‚ ОН : 9?Н ОН 
Ну+аз > а “+аю/2Н. Чт 2 а) На-аз - -® 
9?Н 
;2 4 _ 
На —— д=2 = 0. 


Функция Й/”1(0 определяется из условия И/(1, й, 6) =0: 


И/ (2) =-2) Н, (и,1,1)+ УТ (Аим, (Пк о," 


п=0 = 
Для качественной оценки индуцируемого электрического импульса на его верхней лицевой поверхности 
необходимо сформировать два электрода с радиусом раздела А и подключением их к измерительному прибору. В 
этом случае потенциал фо(х, #), индуцируемый на двух эквипотенциальных поверхностях, представляется в виде: 
фо (и,г) = фи (№)Н (®-к)+ Фе» (КН ("- К), (29) 
где НС...) — единичная функция Хэвисайда. 
Подстановка (29) в (5) позволяет определить выражения для определения потенциалов фо(, фо2(0: 


| Р..отае = | Ра оте = 0. (30) 
0 у 
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В результате решения (30) функции фо(®), фо2(2) определяются следующим образом: 


4и (= 0я[о, (0+0, (+0, $ (д =0я| 0,0 - о, (90.0 


где 


зы Н(п.олуь Чи ЭН5"(т,2и) ав ОНУ (пт) 


а; Ти 02 2-0 Л 02 |2=0 


о 


ти аки, 0) + аи АК» (№=,2) аз АКз (%»,2) | 
: [2=0 


2 т, 
| а5 Л Пе [2=0 Л (и 


: . С ак (АЕ АК (в, 
о-д-® ре абы ий и | 


5 аз 
2 во 22 в || 


1 


Г. 
0: (#)= аз [© (иг) как, О. (1) = ан [© (иг) пак, 
о 


В 


&* Чо (2) 4/1 (2) < Л (1,в) | 410 а: Ч/о (2) аз 4/: (2) 

От 5 | Е а аз Е № +2Ю 2 иги (/,) ое )* ЕЯ а | 
(-)Г 4%() _ 40 ел.) [аю „(ан 4% (2) в (2) 

бы 2 в. 4 [2=0 ь 42 [2=0 - 2 В. д А 4 [2=0 У 42 [5=0 | 


В этом случае разность потенциалов 7(#) определяется равенством: 
Г (1) =Фи (-$%2 (1). (31) 
Результаты исследования. Численные результаты представлены для пластины, изготовленной из 
пьезокерамики состава РИТ-4 [4, 11, 16]: 


{Си, Св, Св, Сзз, С} =13,9;7,78;7,43;1 1,5;2,26} х10'° Па, {Еп,Езз}={6,46;5,62} 10? Ф/м, 
{е15,езл,езз}={12,7; —5,2;15,1} Кл/м*, А=1,6Вт/ (м-К), о, =0,4х10> К", 
к=3х10° Дж/(м?-К), ви =8зз = -0,6х10* КлИм? -К), В» =10* с, “=5,6Вт/(м? -К). 


Исследуется следующий случай изменения температуры ‹ (х,,1.) : 


в! (.1.)= иены „ен |, 


где Тик, (мах — максимальное значение температуры и соответствующее ему время ты пах = 100°С, Ту =20 ие 


На рис.2 приведены графики, отражающие в различные моменты времени мы = 0, 16) изменение 


температуры ®©^(0, 2, 2) в толщине пластины (Б = 14103 м, В" = 1х103 м). 
По результатам расчета наблюдается, что вследствие высокого коэффициента теплопроводности и небольшой 
толщины пьезокерамической пластины, установившийся температурный режим формируется достаточно быстро 


п =10 с) при достижении ©*“(0, 2, 4) на нижней лицевой поверхности (2 = #) 78 °С (рис. 2). 


На рис. 3 показано изменение ©* (0, — ь |] по времени [№ Вот с) с учетом (представлено сплошной 


линией) и без учета (представлено пунктирной линий, В=0) релаксации теплового потока 
(6 = 14х10 5 м, В" = 1х10? м). Необходимо подчеркнуть, что использование гиперболического уравнения 
теплопроводности Лорда-Шульмана необходимо только при исследовании пьезокерамической конструкции 


микро-размеров при очень быстром изменении @э (х,,1,). 
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09° (0, 2,1), °С 


60 


40 


20 


0 0,02 0,04 0,06 2 


Рис. 2. Графики ©” (0,2,#) — 2:1 —=113;2—#=10.%;3—#=1008 


в 
. 
ы 
, 


2х10— 4х10-4 6х10—4 


Рис. 3. Графики 9" (0. ие -(. № 2 


сплошная линия — В=10° (с); пунктирная линия — В=0 


Численные результаты определения функции ©“(, 2, 5) показывают, что при проведении исследования 
конструкции из пьезокерамического материала возможно пренебречь влиянием скорости изменения объема тела 
и напряженности на температурное поле, т.е. использовать в расчетах только уравнение теплопроводности. 

На рис. 4 представлен график перемещений И/\(0, 2, 6) по времени Ё& а на рис. 5 зависимости изменения 
о компоненты нормальных напряжений о»(х, 2,1) по координате г в различные моменты времени: 


1-1=1 


тах> 


— = 10 ах (. и с) ‚ сплошная линия — 2 = 0, пунктирная линия — 2=й. 


Й *(0,2,1)х107,м 


0 2х10* 4х10-4 6х10— 1 


Рис. 4. Графики И” (0, 2, 0) -1 
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сиу, 2, )х10-7, Па 
—1 
2 


0,2 0,4 0,6 0,8 
Рис. 5. Графики ои(у, 2, 8) —г: 1 — {= ах; 2 — = 10ах 


= 


Следует отметить, что в процессе прогрева пластина изгибается с увеличением ее толщины, за счет 
закрепления образуются сжимающие нормальные напряжений с„(г,2,#) во всех точках. В случае полного 
прогрева конструкции (1 = 101„х) величина нормальных напряжений с», 2, 101х) по высоте сечения практически 
совпадают (рис. 5, график 2, сплошная и пунктирная линии). При этом о»(х, 0, 9 остается постоянной на всем 
промежутке времени {> [ах (рис. 5, сплошная линия), а на нижней плоскости в начальный момент времени 
ог, й, й существенно меньше (рис. 5, график 1, пунктирная линия). 

Для качественной оценки индуцируемого электрического импульса в виде разности потенциалов ИКД (31) на 
верхней лицевой поверхности рассматриваемого элемента необходимо сформировать два электрода с радиусом 
раздела А = 0,7 и подключением их к измерительному прибору (рис. 6, сплошная линия). При этом определение 
ИО путем подключения верхней и нижней (заземленной) сплошных электродированных поверхностей пластины 
к вольтметру (рис. 6, пунктирная линия) неэффективно. 


И0х10? 


0,5 


2х10* 4х104 6х10* 


> 
=> 


Рис. 6. Графики /(0 -—1 


Обсуждение и заключение. Разработанное замкнутое решение связанной осесимметричной задачи 
термоэлектроупругости для круглой пластины, выполненной из пьезокерамического материала, является более 
точным по сравнению с решением, которое было разработано при решении задач в несвязанной постановке. 
Связано это с тем, что полученные расчетные соотношения позволяют определить, как влияет на напряженно- 
деформированное состояние и электрическое поле рассматриваемого элемента нестационарное температурное 
поле, что позволяет описывать с большей точностью поведение круглой пьезокерамической пластины под 
действием тепловой и электрической нагрузок. Кроме того, появляется возможность научно обосновать размеры 
двух несвязанных между собой электродов, позволяющих наиболее эффективно измерить индуцируемый 
электрический импульс. 
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